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Abstract 

We complement a recent exact study by L. Samaj on the properties of a guest charge Q immersed 
in a two-dimensional electrolyte with charges +1/ — 1. In particular, we are interested in the be- 
havior of the density profiles and electric potential created by the charge and the electrolyte, and in 
the determination of the renormalized charge which is obtained from the long-distance asymptotics 
of the electric potential. In Samaj 's previous work, exact results for arbitrary coulombic coupling 
(3 were obtained for a system where all the charges are points, provided j3Q < 2 and (3 < 2. Here, 
we first focus on the mean field situation which we believe describes correctly the limit /3 — > but 
(3Q large. In this limit we can study the case when the guest charge is a hard disk and its charge 
is above the collapse value (3Q > 2. We compare our results for the renormalized charge with the 
exact predictions and we test on a solid ground some conjectures of the previous study. Our study 
shows that the exact formulas obtained by Samaj for the renormalized charge are not valid for 
(5Q > 2, contrary to a hypothesis put forward by Samaj. We also determine the short-distance 
asymptotics of the density profiles of the coions and counterions near the guest charge, for arbi- 
trary coulombic coupling. We show that the coion density profile exhibit a change of behavior if 
the guest charge becomes large enough {(3Q > 2 — (3). This is interpreted as a first step of the 
counterion condensation (for large coulombic coupling), the second step taking place at the usual 
Manning-Oosawa threshold f3Q = 2. 

Keywords: Coulomb systems, cylindrical polyelectrolytes, renormalized charge, counterion condensation, 
sine-Gordon model. 
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I. INTRODUCTION 

In a recent paper [lj], Samaj studied the properties of one or two "guest" charges im- 
mersed in a classical (i.e. non-quantum) two-dimensional two-component charge-^mmetric 
electrolyte. Using results from the (l+l)-integrable sine-Gordon model 3, y 
ticular the known expressions for the expectation value of the exponential field 
for the form factors p, 0] of this theory, and the exact solution for the thermodynamics of 
the two-dimensional two-component plasma P| , he was able to determine exactly the excess 
chemical potential of a single "guest" charge immersed in the electrolyte, the long-distance 
behavior of the electric potential created by this guest charge and the long-distance behavior 
of the interaction energy between two guest charges, in the whole regime where the system 
of point charges is stable (i. e. when both the guest charges and the internal charges of the 
electrolytes are point particles). 

An important result from Ref. [lj is for the electric potential created by a single guest 
charge Q immersed in the electrolyte. This potential ip(r) has a long-distance behavior, as 
the distance r — > oo, similar to the screened potential predicted by Debye-Hiickel theory, 

V(r) ~ Q Ieii K (m ir ) (1.1) 

where mi is the inverse screening length (it is also the mass of the lightest breather of 



the sine-Gordon model), and it is given in terms of the inverse Debye length k = y/2ir(3n 
in equation (4.15) of Ref. We shall use the same notations as in Ref. [jj: /3 is the 
Coulombic coupling, the electrolytes charges are +1/ — 1 and n is the density. In (jl.ljl . K 
is the modified Bessel function of order 0. However Q ren is not the charge Q of the guest 
charge (as it would be in Debye-Hiickel theory), but it is known as the renormalized charge. 
Samaj found the following expression for the renormalized charge (equations (5.8) and (5.9) 
of Ref. [l|) 

rwj8/(4-0) tdt ' 

TT/3Q 

(4-/3)sin(^) 

The concept of renormalized charge is very important in colloidal science 1C , 



2 exp 



JO tt sin t , 

Qren = L - 7 SU1 . (1.2) 
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16^ . thus Samaj result is of extreme importance for colloidal science, in particular for the 
study of cylindrical polyelectrolytes, which can be reduced to a two-dimensional problem. 

Samaj claims that the rigorous validity of his result ()1.2j) is for /3\Q\ < 2, which is the 
regime where the system of point particles is stable. However, he gives some arguments to 



support a conjecture he called "regularization hypothesis". This conjecture states that the 
validity of (jl.2j) goes actually beyond (3\Q\ — 2. In the case @\Q\ > 2, the regularization 
hypothesis says that equation (|1.2j) gives the value of the renormalized charge for a guest 
particle of charge Q and radius a in the limit mia <C 1. 

In this article, we present some indications that suggest that the regularization hypothesis 
is not valid. These indications come from the small-coupling limit /3 — *> but when f3\Q\ 
can be arbitrary large. This will be explained in Section ITT1 

In Section IIH[ follows a discussion on the short- distance behavior of the density profiles, 
near the guest charge. In particular, we show that the coion density have a change of 
behavior when (3\Q\ = 2 — j3, which can be interpreted as a "precursor" of the counter ion 
condensation. 

II. THE MEAN FIELD LIMIT: POISSON BOLTZMANN EQUATION 

A. The case a = and (3\Q\ < 2 

Let ip(r) be the electric potential at a distance r from a single guest charge Q immersed in 
the electrolyte. The guest charge is an impenetrable disk of radius a with its charge spread 
over its perimeter. We shall use the dimensionless potential y(f) = /3if>(r) with f = nr. 

For a three dimensional electrolyte in the presence of an arbitrary external charge dis- 
tribution, it is rigorously proved in Ref. [l7[ that, in the limit (3 — > 0, the density and 
correlation functions of the electrolyte are given by the ones of an ideal gas in the presence 
of the external field y(f) which is the solution of the nonlinear Poisson-Boltzmann equation 
with the external source charge. Based on this evidence, we conjecture that this is also 
valid for our problem, although in our case we consider a two-dimensional system, and in 
the case a^Owe include a hard-core interaction between the external guest charge and the 
electrolyte (which is not considered in the proof of Ref. [17]). Thus, assuming the validity 
of this hypothesis, the mean field electric potential y(f) for our problem is the solution of 



A f y(r) = sinh(y(r)) 



r > a 



(2.1) 



A f y(f) = 



r < a 



(2.2) 
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satisfying the boundary conditions 



lim y(f) = (2.3) 

r— >oo 

li m f*V}p_ = _ (3Q / 2 A s 

f^a+ r 

lim y(f) = lim y(f) (2-5) 

where a = na is the guest particle radius in units of 1/k and the charge Q of the guest 
particle is supposed to be uniformly spread over its perimeter. 

This is usually called the mean field approximation. Let us remind the reader that the 
mean field approximation corresponds to the classical treatment of the sine-Gordon model: 
Poisson-Boltzmann equation is the stationary action equation of the sine-Gordon model. Let 
us also clarify, that the nonlinear Poisson-Boltzmann theory is correct in the limit (3 — > 
when it is used to describe the density profiles of the elect roly te created by an external 
charge distribution, under the conditions considered in Ref. |l7[. On the other hand, the 
linear Debye-Hiickel theory should be used to describe the distribution functions of the 



internal charges of the bulk electrolyte in the limit (3 — > 



18] . 



The two-dimensional Poisson-Boltzmann equation (j2.1j) has been solved exactly 



20 



211 ]. and in particular the connexion problem has been studied extensively. This is the 



problem to relate the long-distance behavior of y(f) with its short- distance behavior 



22, 



231 ] . This connexion problem is essential for the determination of the renormalized charge 



Let us recall some of the results from 



21, 



relevant for our discussion. For 



a point guest particle, a = 0, and fl\Q\ < 2, the electric potential has the long-distance 
behavior 

y(r) ~ AXK (f) , f -> oo (2.6) 

and the short- distance behavior 

y n (f) = -2A\nf + 2\nB + o(l) , f -> . (2.7) 



The solution of the connexion problem 
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22j states that the constant A intervening in 



the long-distance behavior is related to the constants A and B of the short-distance behavior 
by 

2 

A = - arcsin(7rA) (2.8) 

7T 
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and 



B 




(2.9) 



where T() is the Gamma function. In relation to the physical problem, one immediately 
recognizes that 2A = (3Q is the bare charge of the guest particle and 4A = /3Q Ten is the 
renormalized charge. Thus, the renormalized charge is given by 



Notice that (j2.10j) corresponds to the limit (3 — > while keeping f3Q arbitrary large (with 
(3\Q\ < 2) of equation (jl.2j) . Indeed the mean field approximation for the study of single 
guest charge Q immersed in the electrolyte is asymptotically correct in the low coupling 
limit (3 — > and (3Q arbitrary. 

B. The case a ^ and (3\Q\ > 2 

If /3 1 Q | > 2 it is mandatory to consider that the guest particle is a disk with impenetrable 
radius a ^ 0, otherwise the charges of opposite sign of the electrolyte will collapse with 
the guest charge. Although it is not possible (yet) to obtain exact results when a ^ for 
arbitrary values of (3, in the limit (3 — > we can obtain some results, under the hypothesis 
that the mean field approach is correct in that limit when o ^ 0. Thus, when a ^ 0, and 
(3 — > and (3\Q\ > 2, we assume that the electric potential is given again by the mean field 
theory: it is the solution of equations (|2.1j) and ()2.2j) and the boundary conditions ()2.Hj) . 
()2.4|) and (|2.5jl . We will consider the case na = a -C 1 but a ^ 0. 

To formally solve this problem, let us introduce yo{r) the solution of Poisson-Boltzmann 
equation (|2.1|) in the whole space, for f e M + , and satisfying the boundary condition ()2.6|) 
when f — ► oo. So far, A in equation ()2.6|) can be seen as an integration constant. Then, the 
electric potential is given by 



I yo(f) for f > a 

Enforcing the boundary condition ()2.4j) at r — a determines the integration constant A. 

To determine A in the case a< lwe need to know the short- distance asymptotics of yo{r). 
Without loss of generality (because the electrolyte is charge symmetric), let us suppose that 




(2.10) 




for f < a 



(2.11) 
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Q > 0. For (3Q < 2 + ©(l/llnal) 1 , the short-distance asymptotics are the same as in the 
previous section, given in equation ()2.7|) . 

If (3Q is large enough (larger than 2) then A > l/V, and then the short-distance behavior 



of yo{r) changes drastically. It is now given by I9I |20| 



Voir) 



-2 In f — ) -2 In 



1 1 

sin I 2/i In — h 20(/i) 



+ 0(f 4 



f -> 



with 



4>{p) = arg(r(l - iji)) 



(2.12) 



(2.13) 



Let us define 

<p{r,n) = 2n\n~ + 24>(ii) (2.14) 

which is the argument of the sine function in ()2.12|) . The constant li appearing in the short- 
distance behavior of yo(f ) is related to A (see equation ()2.16|) below) and can be determined 
by the boundary condition ()2.4|) at r = a. It is the solution of the transcendental equation 



(2.15) 



f3Q = 2 + All cot [2/i In - + 2<j>(ji)\ . 



The solution of the connexion problem gives the following relation between A and fx |19l.l2 



/i = cosh 1 (7rA) 



TT 



(2.16) 



Notice that we choose here /1 < 0. If (3Q > 2 the argument of the cot function in equa- 
tion ()2.15|) is in the range [7r/2,7r]. Equation ()2.16|) allows us to obtain the renormalized 
charge Q ren . 

For a ^ 0, the definition of the renormalized charge comes from the long-distance behavior 
of the electric potential and its comparison with the solution from the linear Debye-Hiickel 
theory 



y[r) 



PQr 



Kn(f) 



T — > OO 



(2.17) 



aKi(a) 

where K\ is the modified Bessel function of order 1. Comparing with ()2.6|) we have j3Q ven = 
4aKi(a)\. Notice the additional factor aKi(a) in the renormalized charge. However this 
factor is not really important since if a 1, aKi{a) ~ 1. We shall see that it is the behavior 



1 When a^O there is a small (negative) correction of order 1/| In a| to the critical value (3Q = 2, for details 
see 0- 
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of A which will be very different from the case when a = and (3Q < 2. Once /x has been 
determined from equation (|2.15jh using equation ()2.16|) we can find the renormalized charge 



PQren = COsh(7T/i) (2.18) 

TT 

Let us comment a few points on the small-f behavior of the potential y(f) and of yo(r). 
From equation (|2.12|) . one can distinguish three special regions. Notice that (p(r,n) is a 
decreasing function of f since \i < 0. The first region is for f = up to a value r* such that 

<p(r*,fj,)=ir. (2.19) 

In this region, the formal solution yo(f) of Poisson-Boltzmann equation has no physical 
meaning, since (p(f,fi) decreases from +oo to tt, and then the argument of the logarithm of 
the second term of equation ()2.12j) oscillates around zero, changing of sign, thus y(f) is not 
always real but can become complex. However this region is inside the guest charge (r* < a) 
and in this region the electric potential is a constant: y(f) = yo{a). 
The second region is for f e [a, tm\ where is given by 

(p(r M ,n)=ic/2. (2.20) 

In this region 7r/2 < ip(f,[i) < tt. If f is close to a (and a is close to r*), the second 
term of y{f) in equation ()2.12|) can be very large because (p(f,fi) is close to tt. Then, as 
f increases, ip(r,fi) decreases from tt down to tt/2, and then, the second term of y(f) in 
equation ()2.12j) decreases very fast. Since the counterion density is proportional to e y ^ r \ 
this indicates that close to the guest charge there is a very large counterion density, which 
decreases quite fast as the distance f increases. This is a manifestation of a widely known 
phenomenon in the theory of cylindrical polyelectrolytes, known as the Manning-Oosawa 
counterion condensation |2jj |26( . The layer of "condensed" counterions extends from r = a 
to f = r M . 

When f = tm the second term of equation ()2.12j) vanishes. At this point the total 
charge Q gue st+condens of the guest particle plus the condensed counterion layer is such that 
/^Qguost+condons = 2, as it can be seen from the first term of equation (I2.12fl . We recover here 
another characteristic of the Manning-Oosawa counterion condensation 25|,|26|,|27|: the layer 



of condensed counterions reduces the bare charge of the guest charge to /3<3 gU est+condens = 2. 
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Above t Mi we enter a third region, outside the condensed layer of counterions, where f starts 
to become large enough such that the small-f behavior (j2.12j) is no longer valid. 

From a physical argument we can now see that as f3Q becomes larger than 2, the renor- 
malized charge is not given anymore by equation (J2.1(Jj) (which is valid only for a = and 
j3Q < 2). Indeed, from the discussion above, we know that in the region r* > r M , just outside 
the counterion condensed layer, the guest charge "dressed" with the condensed counterions 
can be seen as an object with a charge Q gU est+condens = 2//5. Thus its renormalized charge 
will be close to the prediction of equation ()2.10|) for f3Q = 2, that is [3Q ren is close to 4/tt. 

The bare charge could be arbitrary large (with (3Q > 2), but the renormalized charge will 
remain close to 4/tt, because of the Manning-Oosawa counterion condensation phenomenon. 
In particular, the renormalized will not oscillate and become eventually negative as predicted 
by (j2.10j) for (3Q > 2, if the regularization hypothesis was valid. 

This argument can also be justified from a more rigorous point of view. If a is small 
enough, the solution of equation (|2.15J) for /i is small (of order 1/| ln(d/8)|). The renormalized 
charge is given by the correct formula 1)2.18)1 when f3Q > 2. If jj, <C 1 and a <C 1, we see 
from (j2.18j) that f3Q rcn is close to 4/tt. It is actually slightly larger than 4/7T. 

This is verified numerically in figure ^ where we computed the renormalized charge 
from a numerical resolution of Poisson-Boltzmann equation, using the method described 
in Ref. 



28| . We confirm numerically that the renormalized charge, for (3Q > 2, is slightly 
above 4/tc (for a < 1). This can be contrasted with the prediction of the regularization 
hypothesis, shown in dashed line in figure HI where Q mn is expected to decay when [3Q > 2 
and even vanish and change its sign at f3Q = 4. The numerical results show that this is not 
true: Q ren is always an increasing function of Q, and eventually saturates to a finite value 
for large values of Q. 



This saturation phenomenon [2jj of the renormalized charge is quite usual in the nonlinear 
Poisson-Boltzmann approach to the problem. When the saturation phenomenon occurs we 
also have a = r*. Indeed, if a = r*, by the definition (|2.19j) of r* we have (p(a,fi) = n and 
one can verify that in equation (|2.15j) . Q — > +oo. Solving 

(p(a,n)=ir (2.21) 

for /i and replacing in ()2.18j) . allow us to obtain the saturation value of the renormal- 
ized charge. Equation ()2.21|) is a transcendental equation, but since fi is small, of order 
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tea =10 

a=0,pQ<2 

exact result 
__ a=0,pQ>2 

regularization hypothesis 



2 

PQ 



FIG. 1: The renormalized charge Q ren as a function of the bare charge Q, in the mean-field limit 
P — ► 0, for various values of the radius a of the guest charge. For a = 0, the exact result Q2.1U[) is 
shown in full thick line for BQ < 2. In dashed line, the extension of (I2.1UI) for 3Q > 2 is shown: 

n 

this is the prediction of the regularization hypothesis from Ref. j|. The symbols correspond to 
values a > 0, obtained from a numerical resolution of Poisson-Boltzmann equation. 

1/| ln(a/8) |, if a -C 1, it can be solved in an expansion of powers of 1/| ln(a/8)|. For ex- 

I I 

ample, up to order 4 in l/|ln(a/8)|, we find the renormalized charge at saturation [23j] 



(3Q 



sat 
ren 



AaKxia) 



7T 



cosh • 



7T 



+ 0(|lna|- 5 ) 



(2.22) 



2(ln| + 7 ) 
where 7 ~ 0.5772 is the Euler constant. 

In conclusion to this part, we notice the failure of the regularization hypothesis in the 
limit (3 — > with f3Q > 2 and a C 1. The renormalized charge is not given by equation ()1.2|) 
in that limit as the regularization hypothesis claims. 



III. SHORT-DISTANCE BEHAVIOR OF THE DENSITY PROFILES: A "PRE- 
CURSOR" OF THE COUNTERION CONDENSATION AT (3Q = 2 - (3 

An important factor, which is responsible of the failure of the regularization hypothesis 
exposed in the previous section, is the change of behavior of the electric potential at short 
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distances when f3Q > 2. In this section we consider the general situation when < (3 < 2 
and we return to the case of point particles a = 0. We study the short- distance behavior of 
the density profiles and show that there is a change of behavior in the asymptotic expansion 
at short distances of the coion density profiles when f3\Q\ = 2 — (3. 

The short- distance behavior of the density profiles, in the presence of the guest charge, 



can be obtained by adapting an argument presented in Refs. [29|, |30J for the correlation 
functions. Let suppose, without loss of generality, that Q > 0. From the general principles 
of statistical mechanics we know that the short- distance behavior of the density profiles, 
near the guest charge at the origin, is dominated by the Boltzmann factor of the Coulomb 
potential e~^ qVc< ^\ with v c (r) = — lnr the Coulomb potential. We have 

n q {r) ~ n q c Q y Qq , r^O (3.1) 

with q = ±1, n q is the bulk density of charges q, and the constant cg q is related to the excess 
chemical potentials (fi q xc , /ig c , /^f+q) of the charges q, Q and Q + q, which can in turn be 
expressed as expectation values of exponentials of the sine-Gordon field 0, 

( p ib{Q+q)4>\ 

c Qq = eM~m+ q ~ + -/'D] = {elbQ , ){e Z) ( 3 - 2 ) 
with the same conventions as in Q for the normalization of the sine-Gordon field and 
b 2 = /3/4. In the case q = +1 (coions, same sign as Q) this is valid provided f3Q is small 
enough, as we will explain below. 

Let be the grand canonical partition function of the system composed by the elec- 
trolyte and the guest charge Q fixed at the origin, with fugacities z + and ^_ for positive and 
negative particles respectively. The partition function is well defined for point particles if 
f3 < 2 and (3\Q\ < 2. The density of particles of charge q is 

n q (r)=z q r^^^ (3.3) 

where E[Q;q, r] is the partition function of the electrolyte in the electric field created by a 
guest charge Q fixed at the origin and a charge q fixed at r. 

When r — > 0, E[Q; q, r] has a finite limit provided that (3\Q + q\ < 2, since E[Q; q, 0] = 
S[Q + q] is the partition function of a system composed by the electrolyte and a guest charge 
Q + q at the origin. Under this condition we can affirm that the short-distance behavior (|3.1|) 
for the density profile is valid. This can also be seen from ()3.2|) : the expectation value 
(e ib( - Q+ ^) is finite provided (3\Q + q\ < 2 Q. 
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For Q > and q = 1, the above condition reads f3Q < 2 — {3. When (3Q > 2 — (3, 
equation (j3.1|) is no longer valid. We must return to the general expression ()3.3|) . and study 
the short- distance behavior of E[Q;q, r]. If (3Q > 2 — /3, H[Q; g, r] diverges as r — > 0. Its 
short- distance behavior is dominated by the approach of a charge —1 to the system of the 



charges Q and +1 which are separated by a distance r 



29]. 



This gives 



S[g ; g,r] ~ est x r- /3(Q+1)+2 , r^O (3.4) 

with est some constant, which will not be needed for our analysis, but that can eventually 
be evaluated j^ll [3- Thus, the short-distance behavior of the coion density profile, for 

(3Q > 2 - /3, is 

n + (r) ~ est x r 2 ~ p , r — >• (3.5) 

Notice that, on the other hand, the counterion density profile n_(r) behaves always as 
predicted by equation (|3.1|) . since the corresponding E[Q; — l,r] has always a finite limit at 
r = 0, provided /3 < 2 and < 2. 

The mean interaction potential w +i Q(r) of a coion (charge +1) with the guest charge Q 
and its polarization cloud, is defined by n + {r) = n + e~ l3w+ - Q ^ r \ From (|3.5j) . we deduce that 
its short- distance behavior is 

f^Qlnr + 0(1), if/3Q<2-/3 
/3w +>0 (r) = { , r ^ (3.6) 

[(2-/3)lnr + 0(1), if (3Q > 2 - /3 

Notice that, in particular, the short-distance leading behavior of w + ^q{t) is independent of 
Q when f3Q > 2 — j3. This situation can be interpreted as a "precursor" of the Manning- 
Oosawa counterion condensation. When [3Q increases above 2-/3, the counterion cloud 
near the guest charge reduces its bare charge so that the coions "see" a "dressed" object 
with charge 2/(3 — 1, which is independent of Q. 

On the other hand, the counterion mean interaction potential has always the behavior 
(provided (3 < 2 and (3Q < 2) 

(3w_ Q (r) ~/?Qlnr + 0(1), r^O. (3.7) 

The counterions continue to "see" the bare guest charge Q even if (3Q > 2 — (3. 

When (3Q > 2, we arrive at the collapse of the charge Q with the counterions. In this 
case we need to consider that the guest charge is an impenetrable disk with radius o ^ 0. 
(3Q = 2 corresponds to the well-known Manning threshold for counterion condensation. 
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We would like to stress that the above analysis is valid for large coup ling < (3 < 2. 
The usual presentation of the counterion condensation phenomenon j^, I2J is done in a 
small-coupling approximation (3 — > 0. Notice that in this case both limits (3Q < 2 — (3 and 
(3Q < 2 coincide. Here we have put in evidence that, at large coulombic coupling f3, the 
counterion condensation might actually take place in two steps: first when (3Q = 2-/3, 
where the coion density changes its short- distance behavior, but not the counterion density, 
and a second step, at the usual threshold (3Q = 2. 

Let us conclude this section with a conjecture, suggested by Samaj 33^ The mean field 
analysis of the previous section shows that the formula (jl.2j) for the renormalized charge 
is not valid beyond (3Q = 2 in the limit (3 — ■> 0. For arbitrary values of (3, the validity 
of (JH2J) might be beyond f3Q = 2. In the case f3 0, the failure of (JU2D when (3Q > 2 
is accompanied by a change of behavior in the short distance asymptotics of the electric 
potential. 

For arbitrary (3, when (3Q > 2 it is necessary to introduce a hard core for the guest 
particle. The coion density will certainly change its short distance behavior at (3Q = 2. On 
the other hand there are indications that the counterion density will still behave at short 
distances as r _/3( ^ up to (3Q < 2+ (3. This is because, as previously noted, 5[Q; —1, r] remains 
finite up to (3Q < 2 + (3, for small r. Since the counterion density determines the dominant 
behavior at short distances of the electric potential, the electric potential might actually 
not change its short distance behavior until (3Q > 2 + (3. If this is true, the regularization 

n, 

hypothesis put forward by Samaj, might be valid up to (3Q < 2 + (3 33]. Interestingly, 
this leaves open the possibility of charge inversion (i.e. the renormalized charge becomes 
negative) if (3 > 1, when (3Q > 4 — (3. 



IV. CONCLUSION 

As a complement to Ref. [l], we have studied the low-coupling, mean field situation, 
/3 — > 0, but (3\Q\ arbitrary, in order to determine the behavior of the renormalized charge 
when > 2 of a guest charge Q immersed in an electrolyte. We have shown that, at least 
in this mean field situation, the regularization hypothesis put forward by Samaj in Ref. Q] 
is not valid: the formula (jl.2)l for the renormalized charge is not valid when (3\Q\ > 2. 

We have also studied the short-distance asymptotics of the density profiles. The coion 
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density profile exhibits a change of behavior if the guest charge Q is large, (3Q > 2 — /3, as 
shown in equation (|3.6j) . Colloquially speaking, it is like if the coions "see" a "dressed" charge 
2/(5 — 1 instead of Q, when (3Q > 2 — /?. We interpret this as a first step in the Manning- 
Oosawa counterion condensation when the coulombic coupling f3 is large, the second step 
taking place when j3Q = 2 as it is usually explained in the literature |25j, |26J] for the small 
coupling P — > situation. 
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